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Cosmic rays are the most energetic particles arriv- 
ing at earth. Although most of them are thought to 
be accelerated by supernova remnants, the details of 
the acceleration process and its efficiency are not well 
determined. Here we show that the pressure induced 
by cosmic rays exceeds the thermal pressure behind 
the northeast shock of the supernova remnant RCW 
86, where the X-ray emission is dominated by syn- 
chrotron radiation from ultra-relativistic electrons. 
We determined the cosmic-ray content from the ther- 
mal Doppler broadening measured with optical spec- 
troscopy, combined with a proper-motion study in X- 
rays. The measured post-shock proton temperature 
in combination with the shock velocity does not agree 
with standard shock heating, implying that >50% of 
the post-shock pressure is produced by cosmic rays. 

The main candidates for accelerating cosmic rays up 
to at least 10^^ eV are shell-type supernova remnants 
(SNRs), which are the hot, expanding plasma shells, 
caused by exploded stars (supernovas). In order to main- 
tain the cosmic -ray energy density in the Galaxy, about 
3 supemovae per century should transform 10 percent 
of their kinetic energy in cosmic-ray energy. Indeed, 
~ 10^** eV electrons have been detected at forward 
shocks (dl |2ll and possibly at reverse shocks (|3l HI of 



several shell-type remnants by their X-ray synchrotron 
emission, and particles with TeV energies have been de- 
tected in several SNRs by Cherenkov telescopes (|5l[6ll. 

If SNRs transform a substantial amount of their ki- 
netic energy into cosmic rays, this should affect the 
kinematics of the remnant. One imprint of energy losses 
by cosmic rays is a higher compression factor of the 
post-shock plasma(7 ), for which indications have been 
found in both the Tycho SNR and SN 1006 (S ill. An- 
other signature of the energy absorbed by cosmic rays 
is a lower post-shock temperature (fTOl [TTl [121 ITsT l. For 
shocks with conservation of mass, momentum and en- 
ergy, in absence of cosmic rays, the post-shock temper- 
ature (Ti) for species with mass rrii relates to the shock 
velocity (Vg) as 

3 2 

kT, ^ —m,v^ (1) 

in case of no thermal equilibrium (i.e. the several 
atomic species do not have the same temperature), in 
which case protons carry most of the thermal energy. 
In case of fast thermal equilibration, this relation reads 
kT = j^^nipV^ {fx ~ 0.6 for cosmic abundances). In- 
dications for a lower post-shock electron temperature 
have been found in the Magellanic Cloud remnant IE 
0102-72 (14), which may constitute only a minor part 
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Fig. 1: (Left) The eastern rim of RCW 86, as observed in 2007 with Chandra. Red indicates the 0.5-1.0 keV 
band, green the 1.0-1.95 keV band and blue shows the 1.95-6.0 keV band. The northern part has relatively more 
flux in the higher energy bands, which is characteristic for synchrotron emission. (Right) Blue is the broadband 
keV Chandra image, red is the image as observed with the VLT through a narrow Ha filter. The regions (yellow 
and red) indicate where we measured the proper motion. In both panels the location where we took the optical 
spectrum is indicated with a white line. 



of the thermal pressure. Here we derive the post-shock 
proton temperature and the shock velocity of the north- 
east rim of the shell-type SNR RCW 86 based on optical 
and X-ray observations. 

RCW 86 jTSl) was detected in TeV energies by the 
H.E.S.S. telescope (W\ and is probably the remnant of 
the supernova witnessed by Chinese astronomers in 185 
A.D. dTTl [TsT l. It has been suggested that it evolves in 
a stellar-wind blown cavity, where the southwest corner 
has already hit the cavity wall(19). The northeast side 
still expands in a less dense medium and its X-ray spec- 
trum is dominated by synchrotron radiation, which is an 
indication for efficient cosmic-ray acceleration. 

The optical spectrum of the northeast rim of RCW 86 
is dominated by hydrogen lines, with no [Nil] line 
emission(20). The lack of [Nil] indicates that the hy- 
drogen line emission is not a result of strong cooling, 
but results from excitation processes immediately be- 
hind the shock front. The hydrogen lines from these 



shocks consist of two superimposed Gaussian line pro- 
files: one, caused by direct excitation, has the thermal 
width of the interstellar medium (ISM), the other is 
emitted after charge exchange between hot post-shock 
protons and cold incoming neutral hydrogen and hence 
has the thermal width of the post-shock protons. Ha 
emission and efficient cosmic -ray acceleration are likely 
to anti-correlate because incoming neutral species are 
likely to damp plasma waves, which are essential for 
shock acceleration (211) and because cosmic rays escap- 
ing ahead of the shock ionize the surrounding ISM and 
decrease the amount of Ha emission. In RCW 86 the 
Ha emission occurs all along the rim, including, al- 
though with weak emission, the parts coinciding with 
X-ray synchrotron emission, where efficient cosmic-ray 
acceleration is likely to occur (Fig. 1). The only other 
remnant in which Ha emission is seen all along the 
shell, including regions with X-ray synchrotron emis- 
sion, is the SN 1006 SNR(l22ll. 
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The right panel in Figure 1 shows both the Ha and the 
X-ray emission of the northeast rim of RCW 86. The 
Ha emission marks the onset of the X-ray synchrotron 
radiation, which indicates that they are from the same 
physical system. 

In order to measure the proton temperature, we used 
long-slit spectra obtained with the visual and near ultra- 
violet FOcal Reducer and low dispersion Spectrograph 
(FORS2) instrument on the Very Large Telescope (VLT) 
( |23T l. We first imaged the northeast side of RCW 86, 
where the X-ray spectrum is dominated by synchrotron 
emission. Using this image as a guide, we pointed the 
slit at a location where the Ha emission is bright (Fig. 
1, Table SI). 

The spectrum's (Fig. 2) measured full width at a half 
of the maximum (FWHM) is 1100 ± 63 km/s (see SOM 
for further details) corresponding to a cr^, = 467 ± 27 
km/s and implying a post-shock temperature of 2.3 ± 
0.3 keV. 
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Fig. 2: The Ha spectrum, with broad and narrow com- 
ponents (dotted). The best fitting spectrum is overplot- 
ted. The lower panel shows the residuals divided by the 
errors. 



To measure the shock velocity of the northeast rim of 
RCW 86, we observed it with the Chandra X-ray obser- 
vatory in June 2007, and matched it with an observation 
taken in June 2004 dTSl) . To make both observations as 
similar as possible, we used the same observation pa- 
rameters as in 2004 (Table S2). 

We measured the proper motion of the shock at the 
location of the slit of the Ha spectrum by comparing 
the positions of the shock in the two images (see SOM 



for further details). A solid estimate of the proper mo- 
tion is 1.5 ± 0.5" in 3 years time (Fig. 3 and SI), im- 
plying a shock velocity of (6.0 ± 2.0) x 10'^ km/s at 
a distance of 2.5 kpc ([25l l26b . The statistical error on 
the measured expansion is of the order of 0.2". How- 
ever, in the process of calculating the proper motions, 
we found that small details, such as slightly changing 
the angle in which we made the profile, tended to give 
a different proper motion, with a difference larger than 
the 0.2" statistical error we measured. However, in 
none of the measurements, did we find a proper mo- 
tion below 1.0". Because the proper motion is higher 
than expected (18), we verified that it is consistent with 
data taken in 1993 with the Position Sensitive Propor- 
tional Counter (PSPC) on board the ROentgen SATellite 
(ROSAT) compared with the 2007 observation (Fig. 3 
and SI). Although the proper motion, using the nomi- 
nal pointing of the ROSAT PSPC, is statistically highly 
significant, the large pointing error of ROSAT (~4") re- 
sults in a detection of the proper motion at the 2a level. 

Compared to other remnants of a similar age, the 
shock velocity is surprisingly high. Recent models (IZTT i 
predict ^ 5000 kms^^ after 2000 years for SNRs 
evolving in a wind blown bubble (28). This fits with 
the scenario where RCW 86 is evolving in a cavity and 
the southwest corner, which has a slower shock veloc- 
ity (29'; "30^ and a mostly thermal ([3]l X-ray spectrum, 
has already hit the cavity shell. Shock acceleration the- 
ory suggests that only shocks with velocities exceeding 
2000 km s^ ^ emit X-ray synchrotron emission dSTl [TS]) . 
which is also consistent with observations ([32T l. 

An additional uncertainty in the shock velocity is in 
the distance to RCW 86, which is based on converging 
but indirect lines of evidence. RCW 86 was found to be 
in the same direction as an OB association, at a distance 
of 2.5 kpc(33 ). Because high mass stars are often found 
in such associations, the progenitor of RCW 86 may 
well have formed in this one, provided that RCW 86 
is the remnant of an exploded massive star. Other stud- 
ies (I25l,f26^ found a distance of 2.3 and 2.8 kpc respec- 
tively, based on the line-of-sight velocity of ISM swept 
up by the remnant, combined with an observationally 
determined rotation curve of the Galaxy (34). The third 
argument supporting a distance of 2.5 kpc is the molec- 
ular supershell seen in CO emission in the direction of 
RCW 86, whose line-of-sight velocity agrees with that 
of RCW 86(l35Tl. In further calculations, we take the dis- 
tance towards RCW 86 to be 2.5 ± 0.5 kpc, leading to a 
shock velocity of 6000 ± 2800 km/s. 
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Fig. 3: (Left) Steep gradient in the radial profiles for the 2004 and 2007 observations, adaptively binned with the 
Haar method (l24l l. so that each bin has a signal to noise ratio of at least 4, with a maximum bin width of 2". (Right) 
statistics of the proper motion measurement. (Details of the used radial profiles in SOM). 



The relation between shock velocity and measured 
post-shock proton temperature has been extensively 
studied (i20l [30l |36l |32l [ID, including the cross sec- 
tions for excitation and charge exchange as function of 
Vs. Although recent studies show that there can be a sub- 
stantial effect of cosmic rays on the post-shock proton 
spectrum |39]l. up to now, there was no need to include 
cosmic -ray acceleration in the interpretation of the post- 
shock temperature. This is possibly because most of the 
Ha spectra are taken from the brightest rims of SNRs. 
Because Ha emission and efficient cosmic -ray accelera- 
tion are likely to anti-correlate (21), these rims probably 
have low cosmic -ray acceleration efficiency. A possible 
exception is 'knot g' in the Tycho SNR, where indica- 
tions for cosmic-ray acceleration in the form of a pre- 
cursor have been found (l40l l4Tt |42] |. Additionally, for 
some SNRs (" 221 l43l l the distance has been determined 
using the post-shock proton temperature in combination 
with the proper motion, using theoretical models which 
do not take into account energy losses and cosmic-ray 
pressure. This procedure leads to an underestimate of 
the distance if cosmic-ray acceleration is present. Thus, 
unless the distance is accurately determined in an in- 
dependent way, there will be no discrepancy between 
the predicted Vg, based on kT and Eq. [Tjand the actual 
shock velocity. 

The shock velocity of the X-ray synchrotron rim im- 
plies a post-shock temperature of 70 keV (assuming no 
thermal equilibrium), 42 keV (assuming equilibrium), 
whereas the measured post-shock temperature is 2.3 



keV. This measurement is at least a factor 18 less than 
the post-shock temperature estimated from the shock ve- 
locity, which can now be used to constrain current the- 
oretical shock heating models (fT2l [TsT l. Additionally, 
this proton temperature is close to the electron temper- 
ature at the same location(18), implying fast thermal 
equilibration between both species, breaking the trend 
between the shock velocity and the measure of thermal 
equilibrium seen in previous observations! 38l [37]) . 

To translate this discrepancy into the energy and 
pressure in cosmic rays, we followed the approach of 
(11), which is based on standard shock equations for 
plane-parallel, steady-state shocks, modified by addi- 
tional pressure and loss terms [see also ( 1441 |2l l45t l. 
The loss term is defined in terms of the incoming en- 
ergy flux: ecR = FcK/\pavl, Fqk is the amount 
of energy flux in cosmic rays which escapes from the 
system and po is the pre-shock density. The parame- 
ter which indicates the fraction of the pressure induced 
by cosmic rays in the total post-shock pressure is wcR 
(wcR = -Pnt / (^t + ftrT ) , with the pressure in par- 
ticles with a thermal and Pnt with a non-thermal energy 
distribution (i.e. CRs). We plot the modified equations 
(listed in the supporting online material) in Figure 4 and 
indicate the region where the combination of kT and 
Vs of the northern rim of RCW 86 resides for thermal 
equilibrium as well as for no thermal equilibrium. As 
Figure 4 shows, a post-shock temperature and a shock 
velocity do not give a unique solution for lucR and ecR. 
However, the cosmic rays significantly change the shock 
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Fig. 4: The dark blue area is the area allowed according to the modified equations, limited by full thermal equi- 
librium and no thermal equilibrium for the nominal values of kT —2.3 keV, Vs = 6000 km/s. The light blue area 
shows the area of the kT and Vg with all uncertainties taken into account. The thin lines indicate the compression 
ratio (x) of the post-shock plasma. Within the allowed region, A 2 line provides an upper limit on ecR and a 
lower limit on wqr- 



dynamics, because the combination solution is far away 
from wcR — and ecR = (Fig. 4). 

There are two ways to further constrain wqk and ccr- 
First, an additional estimate of the compression ratio 
(x) of the post-shock plasma would exactly determine 
wcR and e^R. For certain SNRs this is done by deter- 
mining the distance between the supernova ejecta and 
the outer shock; a higher compression ratio implies that 
the swept-up ISM forms a thinner shell and hence the 
ejecta will be closer to the shock front (l8l[9]l. However, 
ejecta and swept-up ISM are only distinguished by their 
thermal spectra, which is (almost) absent in the X-ray 
synchrotron dominated rim (~ 15% of the total X-ray 
emission (!T8|). 

An other way is to invoke a dependency of wcr 
on ecR. According to non-linear shock acceleration 
theorv(l46l rnt . ccr/wcr = |(1 - 1/x)^. in which 
X is the compression ratio of the post-shock plasma and 
A = 1, 2 indicates the (wcRj ecr) relation for a cosmic- 
ray spectrum with f{p) oc p^^- , p~^-^ respectively, with 



p the momentum of the cosmic rays. The A = 2 line 
gives an upper limit to the energy losses, since it is 
valid for the most efficient cosmic-ray acceleration by 
cosmic-ray modified shocks (47]). For f{p) oc 
A — ln{pi-aay^/mc)( 12) which can be large and does not 
provide a lower limit to ecR- Taking the A = 2 line 
as an upper limit for ecR , we find a value for wcR of 
> 50%. One remaining question is whether we should 
include the effects of the turbulent magnetic field. The 
average magnetic field pressure in RCW 86 has been es- 
timated to be Pb = BVStt = 2.3 x lO'^dyn cm^^^ 
for a magnetic field of 24 /zG ( flST ). This is an order of 
magnitude below the thermal pressure, which we esti- 
mate to be Pt = nkT 3.7 x 10"^°dyn cm^^, for n = 
O.l jTsT ) and kT = 2.3 keV. In reality, the magnetic field 
pressure may be higher if one takes full account of its 
unknown, turbulent spectrum. 

In summary, our observations show that the post- 
shock temperature of the northeast rim of RCW 86 is 
lower than expected from standard shock relations us- 
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ing the measured shock velocity. The high velocity 
(6000 ± 2800 km/s) of the shock implies a local low 
ISM density, which can be expected in a cavity blown 
by a stellar wind. Cosmic -ray acceleration decreases the 
post-shock proton temperature in RCW 86 by a factor of 
18, implying that > 50% of the post-shock pressure is 
due to cosmic rays. 
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A Supporting online material 

VLT spectra Table [T] lists observational characteris- 
tics of the VLT observation. To reduce the data, we 
first calculated for each pixel the median signal in a 
square of 5x5 pixels around it. Then, the pixels which 
had a higher value than 2x the median were flagged 
as cosmic -ray pixels, and are not used in the following 
calculations. We removed the skylines by fitting a 3"^^ 
order polynomiaQto the spectrum in the spatial direc- 
tion for each wavelength coordinate. Furthermore, we 
calibrated the wavelength using the HeHgCdArNe cali- 
bration spectrum. Then, we added the spectra of the 4 
observations, resulting in the spectrum shown in Figure 
[2] We fitted this spectrum by minimizing the x^, the line 
profile using two Gaussians, convolved with a hat profile 
with the width of the slit (2.5")- The errors on the data 
points are based on the variance of the residuals between 
the data and the best fit model. This method results by 
definition in a x^/d.o.f. ~ 1. The brightness of the 
filaments is 1.0 ± 0.2 x lO^^^erg cm^^ arcsec^^. 

Expansion measurement & statistics Table |2] lists 
the observational parameters of both Chandra obser- 
vations. The data were reduced using the CI AO 
data reduction package, version 4.0 and the calibration 
database CALDB 3.4. We turned off the randomization, 
which is applied during the standard process of gener- 
ating the event list; combined with the dithering of the 
telescope this gives a slightly better angular resolution. 
We checked the registering of the two observations us- 
ing 9 point sources in the field. The error in the align- 
ment of both pointings is well below 0.1", which is the 
value we are using as a systematic error in the proper- 
motion measurement. From the event list, we made ra- 
dial surface brightness profiles with bins in the radial 
direction of 0.25", using photons with energies between 
0.5 and 6.0 keV. We choose the center of the proper mo- 
tion to be a = 14i>42">31.00" and S = -62°29'34.99" 
(J 2000). This is not necessarily the center of the rem- 
nant, it is our estimated center of the curvature of the 
part of the remnant we are interested in. Note that a 
wrongly chosen center can only result in a lower proper 
motion of the shock. We calculated the expansion in 
2 overlapping regions (Figure [TJ. We implemented the 
Poisson statistics as follows: for each bin, we calculated 
the probability that the number of counts in both bins 

' Or a lower order, depending on whether the fit improved by taking 
the higher order. 



were drawn from the same Poisson distribution. We 
measure the proper motion by shifting the normalized 
profiles with respect to each other, with steps of 1 bin 
and for each shift, we calculate the probability that both 
profiles were drawn from the same distribution, using 
the ™d the Poisson maximum likelihood method. 
Additionally, we used the Kolmogorov Smirnov statis- 
tic. For applying the latter, we first made cumulative 
distributions for the profiles and then calculated the Kol- 
mogorov Smirnov statistic for each shift. 

The Kolmogorov-Smimov statistic does not provide 
an error on the parameters. To estimate the errors, we 
simulated 20 radial profiles using the bootstrap method 
( l48l l. added an artificial proper motion (1.0") and 
measured this proper motion using the Kolmogorov- 
Smimov statistic. The best fit proper motion agrees well 
with the input value. We use the standard deviation on 
these best fit proper motions as the la error. In addition, 
we checked the Poison maximum likelihood estimation 
and statistics in a similar way, resulting in consistent 
estimates of the artificial proper motion (1.0" ± 0.2"). 
Furthermore, we use these simuations to determine the 
optimum choice for the range in radius. The most reli- 
able results for the proper motion were obtained if we 
include the whole rim in the radial profile. We vali- 
dated our result using two independently written com- 
puter codes. 

kT and Vg relation Analytical equations of the post- 
shock temperature and shock velocity in the presence of 
cosmic-ray acceleration have been described by several 
authors (I44l [7ll49ll. In this section, we give a summary 
of the equations used to obtain Figure |4] We start out 
with a relation which states the conservation of momen- 
tum over the shock front: 

P2 + P2ul = Po + Poul (2) 

A '0' subscript means pre-shock and a '2' subscript 
means post-shock, p denotes the density, P the pres- 
sure and u the velocity of the gas/plasma in the frame 
of the shock. Now, we use conservation of mass over 
the shock: pqUq = P2U2 and we define the compression 
ratio X = P2/P0: 

P2 = Pa + Poul{l-l/x)- (3) 

We introduce wcR = -Pnt/(-Pt + -Pnt) ^ P2 = 
Pt/(1 — wck)- Additionally, we use Pt ~ nkTi, 
p = nrrii and we assume that the pre-shock pressure 
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Table 1 : Journal of the VLT long slit spectroscopic observations 



Pointing 



a 

(J2000) 



(J2000) 



observation date 



exposure time 



position angle 
(East of North) 



NE 



14:45:15.7 
14:45:15.7 
14:45:15.7 
14:45:15.7 



-62:16:33.2 
-62:16:33.8 
-62:16:34.6 
-62:16:34.0 



05/16/2007 
07/16/2007 
07/18/2007 
07/20/2007 

^Coordinates given are the center of the CCD with which we took the spectrum. The slit, which was not in the 
middle of the CCD, was laid over the filament with coordinates a = 14:45:02.813, S =-62:16:33.05 (J 2000). 



2734 s 
2734 s 
2734 s 
2734 s 



240° 
240° 
240° 
240° 



is small compared to PqUq{1 — 1/x)' 

m = {l-wcK)-il-l/x)mrul (4) 
X 

Note that for a shock without cosmic-ray pressure 
(wcR — 0) and no energy losses and a non-relativisitic 
gas (x — 4), we get equation [T] To derive the compres- 
sion ratio X, we use 7s ; the effective adiabatic index at 
the shock front. This is defined as (i44i) : 

5 + 3wcR 
Is = ^7^— (5) 

3(1 + WCR) 

We now use 7^ in the equation for compression ratio 
including energy losses by cosmic rays, as described in 
(l45l l and define G — |wcR + §■ This gives: 

G+v/G^-(l-ecR)(2G-l) 

X= 1 • (6) 

1 - ecR 
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Table 2: Journal of the Chandra observations 



ObsID 



a 

(J2000) 



S 

(J2000) 



observation start date 



exposure time 



Roll angle 



4611 
7642 



14:45:03.60 
14:45:04.48 



-62:21:05.56 
-62:20:40.53 



06/15/2004 
06/20/2007 



69.1 ks 
71.7 ks 



295.16 
299.02 



Kolmogorov— Smirnov 



Poisson maximum likelihood 



0.03 




0.02 r. 



-0.4 -0.2 0.0 0.2 0.4 0.5 0.! 
shift ["/yr] 



144°- ISO.? 
ROSAT 

other center 



-0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 
shift [Vyr] 



Fig. 5: Plotted are the statistics on the proper motion of RCW 86, for both the Kolmogorov Smirnov statistic and 
the Poisson maximum likelihood estimation. The linestyles denote the same regions as in Figure [3] Since the 
Poisson statistics is a multiplication of the probabilities for each bin, the total probability depends on the number 
of bins in the profile. We corrected for that by subtracting Zn(P,„ax), we multiplied the y-axis with -1, so the best 
fit solution is a minimum in this plot. 
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